SYSTEM AND METHOD FOR PREDICTING 
CHROMOSOMAL REGIONS THAT CONTROL PHENOTYPIC TRAITS 



CROSS REFERENCE TO RELATED APPLICATION 



This application is a continuation-in-part of United States Application no. 
09/737,918, filed December 15, 2000, which is incorporated by reference herein in i 
entirety. 



COMPUTER PROGRAM USTING APPENDIX 



One compact disc that includes a Computer Program Listing Appendix has 
been submitted in duplicate in the present application. The size of the files contained 
in the Computer Program Listing Appendix, their date of creation, their time of 
creation, and their name are found in Table 1 below. In Table 1, each row represents 
a file or directory. If the row represents a directory, the designation "<DIR>" is 
provided in column one. If the row represents a file, the size of the file in bytes is 
provided in column one. Colunms two and three respectively represent the date and 
time of file or directory creation while the fourth column represents the name of the 
file or directory. 

TABLE 1 



Table 1. Contents of the Computer Program Listing Appendix 



Size 


Date of 

Creation 


Time of 
Creation 


File Name 


32,366 
33,987 
11,837 
6,583 
5,141 


12-10-01 
12-10-01 
12-10-01 
12-10-01 
12-10-01 


3.35pm 
3:36pm 
3:36pm 
3:36pm 
3.35pm 


Digidisease.pl 

Display_dev.pm 

Input_output_dev.pm 

Locus_matrix.pl 

Matrix_eenerator.pl 



The Computer Program Listing Appendix disclosed in Table 1 is hereby 
incorporated by reference. 
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Background of the Invention 

Identification of genetic loci that regulate susceptibility to disease has 
promised insight into pathophysiologic mechanisms and the development of novel 
therapies for common human diseases. Family studies clearly demonstrate a heritable 
5 predisposition to many common human diseases such as asthma, autism, 

schizophrenia, multiple sclerosis, systemic lupus erythematosus, and type I and type n 
diabetes mellitus. For a review, see Risch, Nature 405, 847-856, 2000. Over the last 
20 years, causative genetic mutations for a number of highly penetrant, single gene 
(Mendelian) disorders such as cystic fibrosis, Huntington's disease and Duchene 
10 muscular dystrophy have been identified by linkage analysis and positional cloning in 
human populations. These successes have occurred in relatively rare disorders in 
which there is a strong association between the genetic composition of a genome of a 
f^^^ species (genotype) and one or more physical characteristics exhibited by the species 

Q (phenotype). 

15 It was hoped that the same methods could be used to identify genetic variants 

hJ 1 

1*^^ associated with susceptibility to common diseases in the general population. For a 

XI review, see Lander and Schork, Science 265, 2037-2048, 1994. Genetic variants 

associated with susceptibility to subsets of some common diseases such as breast 

|1J cancer (BRCA-1 and -2), colon cancer (FAP and HNPCC), Alzheimer's disease 

Li 

20 (APP) and type n diabetes (MODY-1 , -2, -3) have been identified by these methods, 
which has raised expectations. However, these genetic variants have a very strong 
effect in only a very limited subset of individuals suffering from these diseases 
(Risch, Nature, 405, 847-856, 2000). 

Despite considerable effort, genetic variants accounting for susceptibility to 

25 common, non-Mendelian disorders in the general population have not been identified. 
Since multiple genetic loci are involved, and each individual locus makes a small 
contribution to overall disease susceptibility, it will be quite difficult to identify 
common disease susceptibility loci by applying conventional linkage and positional 
cloning methods to human populations. Mapping of disease susceptibility genes in 

30 human populations has also been hampered by variability in phenotype, genetic 
heterogeneity across populations, and uncontrolled environmental influences. The 
variable reports of linkage between the chromosome lq42 region and systemic lupus 
erythematosus illustrate the difficulties encountered in human genetic studies. One 
group reported strong linkage between the lq42 region (Tsao, J.Clin.Invest, 99, 725- 
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731, 1997) and to microsatellite alleles of a gene (PARP) within that region (Tsao, 
J.Clin.Invest. 103, 1135-1140, 1999). In contrast, no evidence for association with 
the PARP microsatellite marker was noted (Criswell et al., J.Clin.Invest, Jun;105, 
1501-1502, 2000; Delrieu et al.. Arthritis & Rheumatism 42, 2194-2197, 1999); and 
5 minimal (Mucenski, et al., Molecular & Cellular Biology 6, 4236-4243, 1986) or no 
linkage (Lindqvist, et al., Joumal of Autoinmiunity, Mar;14, 169-178, 2000) to the 
lq42 region was found in several other SLE populations analyzed. It is likely that 
additional tools and approaches will be needed to identify genetic factors underlying 
common human diseases. 
10 Analysis of experimental murine genetic models of human disease biology 

should greatly facilitate identification of genetic susceptibility loci for common 
human diseases. Experimental murine models have the following advantages for 
genetic analysis: inbred (homozygous) parental strains are available, controlled 
gl breeding, common environment, controlled experimental intervention, and ready 

^ 15 access to tissue. A large number of murine models of human disease biology have 
been described, and many have been available for a decade or more. Despite this, 
relatively limited progress has been made in identifying genetic susceptibility loci for 
I complex disease using murine models. Genetic analysis of murine models requires 

PJ generation, phenotypic screening and genotyping of a large number of intercross 

20 progeny. Using currently available tools, this is a laborious, expensive and Xxmc- 
consuming process that has greatly limited the rate at which genetic loci can be 
identified in mice, prior to confirmation in humans. For a review, see Nadeau and 
Frankel, Nature Genetics Aug; 25, 381-384, 2000. 

The difficulties encountered in associating phenotypic variations, such as 
25 susceptibility to common diseases, with genetic variations gives rise to a need in the 
art for additional tools for identifying chromosomal regions that are most likely to 
contribute to quantitative traits or phenotypes. In view of this situation, it would be 
highly desirable to provide a technique for associating a phenotype with one or more 
candidate chromosomal regions in the genome of an organism without reliance on 
30 time consuming techniques such as cross breeding experiments or laborious post-PCR 
manipulation. 
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Summary of the Invention 

The present invention provides a system and method for associating a 
phenotype with one or more candidate chromosomal regions in the genome of an 
organism. In the method, phenotypic differences between a pluraHty of strains of the 
5 organism are correlated with variations and/or similarities in the respective genomes 
of the plurality of strains of the organism. The invention relies on the use of a 
genotypic database that includes variations and similarities of representative strains of 
the organism of interest. Representative genotypic databases include, but are not 
limited to, single nucleotide polymorphism databases, microsatellite marker 
10 databases, restriction fragment length polymorphism databases, short tandem repeat 
databases, sequence length polymorphism databases, expression profile databases, 
and DNA methylation databases. 

One embodiment of the present invention provides a method for associating a 
§ phenotype with one or more candidate chromosomal regions in a genome of an 

f :j 15 organism. In this method, a phenotypic data structure that represents a difference in 
one or more phenotypes between different strains of the organism is derived. In its 
I] simplest form, the phenotypic data structure comprises a definition of one or more 
phenotypes exhibited by the organism together with a measure of each of these 
phenotypes. For example, a hypothetical phenotypic data structure for rabbits could 



9 

m 

20 include the phenotypes "tail length" and "hair color" and the respective measure for 



25 



P each of these phenotypes could be "7 centimeters" and "brown. 

A genotypic data structure is established in accordance with one embodiment 
of the present invention. The genotypic data structure is identified by a particular 
locus selected from a plurality of loci present in the genome of the organism. The 
genotypic data structure includes one or more positions within the locus. For each of 
these positions, the genotypic data structure provides information on the extent of a 
variation between different strains of the organism. A hypothetical example of a 
genotypic data structure in accordance with the present invention is a data structure 
for a locus that includes genes A and B. In such an example, the genotypic data 
30 structure includes the positions of genes A and B within the locus as well as some 
measurement related to genes A and B, such as the mRNA expression level that has 
been measured for each of these genes. In this example, the mRNA expression-level 
defines the extent of variation between different strains of the organism. 
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The phenotypic and genotypic data structures are then compared to form a 

correlation value. The process continues with the establishment of another genotypic 

data structure that corresponds to a different loci and the concomitant comparison of 

this genotypic data structure to the phenotypic structure until several of the loci in the 

5 genome of the organism have been tested in this manner. In this way, one or more 

genotypic data structures are identified that form a high correlation value relative to 

all other genotypic data structures that have been compared to the phenotypic data 

structure. Further, the loci in the genome of the organism that correspond to the 

highly correlated genotypic data structures represent one or more candidate 

10 chromosomal regions that may be associated with the phenotype of interest. 

In some embodiments of the present invention, each element in a phenotypic 

data structure represents a variation in the phenotype between a different first and 

I'* second strain of the organism of interest. Such variations may be determined by 

Q measurement of an attribute corresponding to the phenotype in the respective strains 

15 of the organism. Representative phenotypic variations include, for example, eye 

l*-^ color, hair color, and susceptibility to a particular disease. In other embodiments, 

CI 

\| each element in a phenotypic data structure represents a variation in the phenotype 

* between a different first and second cluster of strains of the organism of interest. 

P I In additional embodiments of the present invention, the genotypic data 

20 structure represents a variation of at least one component of a locus between two 
strains of the organism of interest. In other embodiments, each element in the 
genotypic data structure represents a variation of at least one component of the locus 
between a different first cluster of strains of the organism and a different second 
cluster of strains of the organism. In some embodiments, the phenotypic and 
25 genotypic data structures represent a subset of all strains of the organism of interest. 
The present invention contemplates a considerable number of different 
methods for comparing the phenotypic and genotypic data structures. In one 
embodiment the correlation value between the phenotypic data structure and a 
particular genotypic data structure is formed in accordance with the expression: 



a 



30 



S'(p(i) (g(0 -<G^>) 



{ [ r (P(i) - < P >)'] [r (g(i) - < >) '] } 



2t 1 1/2 
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where, 

c(F, G^) is the correlation value; 

p(i) is a value of the i"*" element of the phenotypic data structure; 
g(i) is a value of the i"*" element of the genotypic data structure; 
<P> is a mean value of all elements in the phenotypic data structure; 
<gS is a mean value of all elements in the genotypic data structure; 



and 



5^' =X .where 



N is equal to a number of elements in the genotypic data structure. 
10 Other methods for forming a correlation value between the phenotypic data structure 
and a particular genotypic data structure include but are not limited to regression 
analysis, regression analysis with data transformations, a Pearson correlation, a 
Spearman rank correlation, a regression tree and concomitant data reduction, partial 
least squares, and canonical analysis. 
h'^ .5 In some embodiments of the present invention, statistical methods are used to 

\l identify which of the genotypic data structures that have been compared to a 

l.^ phenotypic data structure are highly correlated. In one such embodiment, a mean 

PJ correlation value that represents a mean of correlation values is computed between the 

phenotypic data structure and a particular genotypic data structure. Further, a 
'^'^ 20 standard deviation of the mean correlation is computed. Genotypic data structures 
having a correlation value that is a number of standard deviations above the mean 
coirelation value are considered to be the data structures that correspond to loci that 
are associated with the genotypic trait. The number of standard deviations that is 
chosen for the cutoff is dynamically chosen so that a specific percentage of the 
25 genome, such as ten percent, is identified as positive. 

Another aspect of the present invention provides a method of determining a 
portion of a genome of an organism that is responsive to a perturbation. In this aspect 
of the present invention, a first phenotypic data structure is produced that represents a 
difference in a first phenotype between different strains of the organism. The first 
30 phenotype is measured for each of the different strains of the organism when each 
different strain is in a first state. Then, a genotypic data structure is established. The 
genotypic data structure corresponds to a locus selected from a plurality of loci within 
the genome of the organism. Further, the genotypic data structure represents a 
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variation, between different strains of the organism, of at least one component of the 
selected locus. The first phenotypic data structure is compared to the genotypic data 
structure to form a correlation value. These establishing and comparing steps are 
repeated for each locus in the plurality of loci. In this way a first set of genotypic data 
5 structures is identified thats form a high correlation value relative to all other 
genotypic data structures evaluated in iterations of the comparing step. 

Then, a second phenotypic data structure is constructed that represents a 
difference in a second phenotype between different strains of the organism. The 
second phenotype is measured for each of the different strains of the organism when 
10 each of the different strains are in a second state that is produced by exposing the 
different strains of the organism to a perturbation. The second phenotypic data 
structure is correlated to the genotypic data structure to form a correlation value. The 
1=* computing and correlating steps are repeated for each locus in the plurality of loci 
§5 thereby identifying a second set of genotypic data structures that forms a high 
!t 15 correlation value relative to all other genotypic data structures that are evaluated 
U during the correlating step. Finally, a dissimilarity in the first set of genotypic data 

structures and the second set of genotypic structures is resolved, thereby determining 
the portion of the genome of the organism that is responsive to the perturbation. 
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O Brief Description of the Drawings 

HG. 1 illustrates a computer system for associating a phenotype with one or more 
candidate chromosomal regions in a genome of an organism in accordance with one 
embodiment of the present invention. 
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HG. 2 illustrates the processing steps for associating a phenotype with one or more 
candidate chromosomal regions in a genome of an organism in accordance with one 
embodiment of the present invention. 

FIG. 3 illustrates a hypothetical representation of the method for computational 
prediction of QTL intervals in accordance with one embodiment of the present 
invention. 
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FIGS. 4A - 4D illustrate the computational prediction of chromosomal regions 
containing genes that determine MHC haplotype (FIG. 4A), lymphoma susceptibility 
(FIG. 4B), airway hyperresponsiveness (FIG. 4C) and retinal ganglion number (FIG. 
4D) in accordance with one embodiment of the present invention. 

FIG. 5 illustrates an analysis of the sensitivity of the computational genome scanning 
method for prediction using ten experimentally verified QTL intervals. A graph of 
the percentage of correct predictions as a function of the amount of genomic sequence 
(percent) contained within the predicted regions is plotted. 



HG. 6 illustrates the comparison of a genotypic database 52 that includes SNP data 
versus a genotypic database that includes microsatellite data in identifying the murine 
chromosomal location for the phenotypic trait of retinal ganglion cell formation, in 
g'l accordance with one embodiment of the present invention. 

15 

m 

h'^ Fig. 7 illustrates a graphical user interface having a toggle that is set to a mode in 

m 

%| which each accession number in a locus L contributes equally to a correspondmg 

genotypic matrix G. 

m 

20 Fig. 8 illustrates a graphical user interface in which a toggle is used to toggle between 
a mode in which each locus position x contributes equally to a corresponding 
genotypic matrix ("by SNP") and a mode in which each accession number contributes 
to a corresponding genotypic matrix. 

25 Fig. 9 illustrates a graphical user interface in which a toggle is provided for switching 
between a weighted mode, in which each computed correlative measure is weighted 
by the number of locus positions x that are represented by a correlative measure, and 
an unweighted mode, in which each computed correlation coefficient is not weighted 
by the number of locus positions x within the respective locus L. 



Fig. 10 illustrates a graphical user interface in which, a user toggle is provided for 
allowing a user to determine the size of the locus L that is used in various 
computations in accordance with one embodiment of the present invention. 
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Fig. 1 1 shows a correlation map where each variation in each locus used to compute 
the correlation map is allowed to fully contribute to the corresponding genotypic 
matrix irrespective of whether multiple variations exists in a single gene. 

5 Fig. 12 shows a correlation map where each gene that includes a variation contributes 
equally to the corresponding genotypic matrix irrespective of the number of variations 
in each gene. 

Like reference numerals refer to corresponding parts throughout the several 
10 views of the drawings. 



C 



Detailed Description of the Invention 

A key aspect of research in genetics is associating sequence variations with 
Q heritable phenotypes. The most common variations are single nucleotide 

ylj 15 polymorphisms (SNPs), which occur approximately once every 100 to 300 bases in a 

genome. Because SNPs are expected to facilitate large-scale association genetics 
|| studies, there has recently been great interest in SNP discovery and detection. The 

f present invention contemplates the use of genotypic databases such as SNP databases 

P.! in order to correlate genetic variances in an organism with one or more phenotypic 

20 variances. As an example, a searchable database of mouse SNPs that contains alleles 
^ for 15 conmion inbred mouse strains and information for performing high throughput, 

inexpensive genotyping assays for each SNP was built. Using pooled DNA samples 
and SNP genotyping assays in the database, a genome scan on phenotypically extreme 
progeny from an experimental intercross was completed. SNP-based genotyping of 
25 pooled samples requires at least twenty-fold fewer assays than genotyping individual 
samples with microsatellite markers, and identified the same linkage regions. 

Although the examples provided herein utilize a genotypic database that 
includes fifteen mouse strains, it will be appreciated that the methods of the present 
invention allow for the use of any number of different types of genetic information. 
30 For example, suitable genotypic databases include databases that have various types 
of gene expression data from platform types such as spotted microarray (microarray), 
high-density oligonucleotide array (HDA), hybridization filter (filter) and serial 
analysis of gene expression (SAGE) data. Another example of a genetic database that 
can be used is a DNA methylation database. For details on a representative DNA 
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methylation database, see Grunau et a/., "MethDB- a public database for DNA 
methylation data," Nucleic Acids Research, in press; or the URL: http://genome.imb- 
jena.de/public.html. 

Gene expression changes often reflect genotypic variation. Therefore, 

5 databases of gene expression among tissues obtained from different individuals 
(mouse strains or humans), can also be utilized by this method. The chromosomal 
position of all human genes is known for human genes, as a result of physical 
mapping or sequencing of the human genome. For gene expression data for mouse or 
other species, the chromosomal location is either known (physical mapping or mouse 

10 genomic sequencing) or can be estimated by syntenic mapping based upon homology 
with human genes. 

To accelerate the process of analyzing experimental genetic models in order 
to identify the genetic causes of complex human disease, the present invention 



C3 provides tools for scanning genotypic databases, such as SNP databases, to predict 

Ul 15 quantitative trait loci (QTL) after phenotypic information obtained from common 
pi Strains of the organism is provided. The computational QTL prediction method is 

M capable of correctly predicting the chromosomal regions that have been previously 

identified by tedious and laborious analysis of experimental intercross populations for 
Pf the multiple traits that are analyzed. Thus, the present invention bypasses the 

20 burdensome requirement for generation and characterization of intercross progeny, 

Q 

enabling QTL regions to be predicted within a millisecond time frame. 

FIG. 1 shows a system 20 for associating a phenotype with one or more 
candidate chromosomal regions in a genome of an organism. 

System 20 preferably includes: 
25 •a central processing unit 22; 

• a main non-volatile storage unit 34, preferably a hard disk drive, for 
storing software and data, the storage unit 34 controlled by disk 
controller 32; 

• a system memory 38, preferably high speed random-access memory 
30 (RAM), for storing system control programs, data, and application 

programs, including programs and data loaded from non-volatile 
storage unit 34; system memory 38 may also include read-only 
memory (ROM); 
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• a user interface 24, including one or more input devices (26, 30) and a 
display 28; 

• a network interface card 36 for connecting to any wired or wireless 
communication network; and 

5 •an internal bus 33 for interconnecting the aforementioned elements of 

the system. 

Operation of system 20 is controlled primarily by operating system 40, which 
is executed by central processing unit 22. Operating system 40 may be stored in 
system memory 38. In a typical implementation, system memory 38 includes: 

• operating system 40; 

• file system 42 for controlling access to the various files and data 
structures used by the present invention; 

• phenotype / genotype processing module 44 for associating a 
phenotype with one or more candidate chromosomal regions in a 
genome of an organism; 

• genotypic database 52 for storing variations in genomic sequences of a 
plurality of strains of an organism; and 

• phenotypic data 60 that includes measured differences in one or 
phenotypic traits associated with the organism. 

In a preferred embodiment, phenotype / genotype processing module 44 
includes: 

• a phenotypic data structure derivation subroutine 46 for deriving a phenotypic 
25 data structure that represents a variation in a phenotype between different 

strains of an organism of interest; 

• a genotypic data structure derivation subroutine 48 for establishing a 
genotypic data structure that corresponds to a locus in the genome of the 
organism of interest; and 

30 •a phenotype / genotype comparison subroutine 50 for comparing the 

phenotypic array to the genotypic array to form a correlation value. 
The operation of these subroutines is described below in connection with the 
description for Fig. 2. 

II 288881.6 
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Genotypic database 52 is any type of genetic database that tracks variations in 
the genome of an organism of interest. Information that is typically represented in 
genotypic database 52 is a collection of loci 54 within the genome of the organism of 
interest. For each locus 54, strains 56 for which genetic variation information is 
available are represented. For each represented strain 56, variation information 58 is 
provided. Variation information 58 is any type of genetic variation information. 
Representative genetic variation information 58 includes, but is not limited to, single 
nucleotide polymorphisms, restriction fragment length polymorphisms, microsatellite 
markers, restriction fragment length polymorphisms, and short tandem repeats. 
Therefore, suitable genotypic databases 52 include, but are not limited to: 



Genetic 
variation type 


Uniform resource location 




http://bioinfo.pal.roche.com/usuka_bioinformatics/cgi- 
Din/msnp/msnp.pi 


SNP 


http://snp.cshl.org/ 


SNP 


http://www.ibc.wustl.edu/SNP/ 


SNP 


http://www-genome.wi.mit.edu/SNP/mouse/ 


SNP 


http://www.ncbi.nlm.nih.gov/SNP/ 


Microsatellite 
markers 


http://www.informatics.jax.org/searches/polymorphism_form.shtml 


Restriction 
fragment 
length 
polymorphisms 


http://www.informatics.jax.org/searches/polymorphism_form.shtmI 


Short tandem 
repeats 


http://www.cidr.jhmi.edu/mouse/nmiset.html 


Sequence 
length 
polymorphisms 


http://mcbio.med.buffalo.edu/mit.html 


DNA 
methylation 
database 


http://genome.imb-jena.de/public.html 



In addition, the genetic variations used by the methods of the present invention 
may involve differences in the expression levels of genes rather than actual identified 
variations in the composition of the genome of the organism of interest. Therefore, 
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genotypic databases 52 within the scope of the present invention include a wide array 
of expression profile databases such as the one found at the URL: 

http://www.ncbi.nlm.nih.gov/geo/ 
It will be appreciated that when the variation tracked by genotypic database 52 is a 
variation in the expression level of a gene rather than a variation in the genome, there 
is no requirement that genomic database 52 be populated with elements such as locus 
54. 

Referring to FIG. 2, the processing steps that are performed in accordance 
with one embodiment of the present invention are illustrated. In processing step 202, 
a phenotypic data structure is derived from phenotypic data 60 (FIG. 1) using 
phenotypic data structure derivation subroutine 46 (FIG. 1). The phenotypic data 
structure tracks measured differences in traits between strains of an organism of 
interest. 

In one embodiment, the phenotypic data structure used is a phenotypic array. 
In this embodiment, the phenotypic array is formed in a stepwise fashion by 
subroutine 46. First, an W x AJ phenotypic distance matrix, P, is established where 
both the ah row and the Jih column are associated with a given strain for which 

quantitative information tj exists for a given trait. 

This matrix is populated with the differences between strains in regard to the 
examined trait as follows: 

Therefore, each element in the matrix corresponds to a distance between strains using 
the quantitative trait as a metric for the space. This matrix has the following 
properties: 

• All of its diagonal elements are zero, because 

p(i6 0^\t^-tj =0 V 

• The matrix is symmetric, because 
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As an example, consider phenotypic information on the lifespan of five mouse strains: 



Strains 


Lifespan (days) 


A /T 

A/J 


til 


A \^X> /T 

/\JvK/J 


ZoZ 


C3H/HeJ 


510 


C57BI76J 


895 


DBA/2J 


568 



1=* 



An exemplary phenotypic distance matrix that tracks the lifespan for these five 
species members has the form: 



p 


A/J 


AKR/J 


C3H/HeJ 


C57BLy6J 


DBA/2J 


A/J 


0 


495 


267 


118 


209 


AKR/J 


495 


0 


228 


613 


286 


C3H/HeJ 


267 


228 


0 


385 


58 


C57BL/6J 


118 


613 


385 


0 


327 


DBA/2J 


209 


286 


58 


327 


0 



10 



Each value in this illustrative phenotypic distance matrix represents the difference in 
life span between the designated members. 

The phenotypic data structure derivation subroutine 46 converts the 
phenotypic matrix to the phenotypic array by taking the non-redundant, non-diagonal 
elements of the matrix and arranging them into a vector P\ 



P^piiaX p(l,3), p(l,N), p(2,3), p(2,4), p(2,N), ...p(N-l,N) 

The vector P obtained for the illustrative distance matrix set forth above is P = (495, 
15 267, 118, 209, 228, 613, 286, 385, 58, 327). The linear format of P facilitates the 

ordered comparison of the phenotype and genotype of respective strains of an 

organism of interest in subsequent computational steps. 

In some embodiments of the present invention, the phenotypic data used by 

phenotypic data structure derivation subroutine 46 (FIG. 1) in processing step 202 
20 (FIG, 2) is entered by hand into system 20 by a computer operator. In other 

embodiments, the phenotypic data is read from a source such as phenotypic data file 



14 



288881.6 



60 (FIG. 1). It will be appreciated that there are no Hmitations on the format of the 
phenotypic data. The phenotypic data can, for example, represent a series of 
measurements for a quantifiable phenotypic trait in a collection of strains of a species. 
Such quantifiable phenotypic traits may include, for example, murine tail length, 
5 hfespan, eye color, size and weight. Alternatively, the phenotypic data can be in a 
binary form that tracks the absence or presence of some phenotypic trait. As an 
example, a "1" may indicate that a particular species of the organism of interest 
possesses a given phenotypic trait and a "0" may indicate that a particular species of 
the organism of interest lacks the phenotypic trait. The phenotypic data structure can 
10 be populated with any form of biological data that is representative of the phenotype 
of the organism of interest. Thus, in some embodiments of the present invention, the 
phenotypic data can be expression data such as mRNA expression data or protein 
1"^ expression level data. In such embodiments, each element in the phenotypic data 

P structure is populated with differences in mRNA or protein expression levels between 

ft 15 Strains of the organism of interest or of cells cultured from the organism of interest. 
h'^ In processing step 204, a particular locus is selected withm the genome of the 

\| organism of interest. Processing step 204 is the first step of a repetitive loop formed 

^ by processing steps 204 through 212 that is repeated for several different loci, or 

fij positions, within the genome of the organism of interest. In some embodiments of the 

tit 20 present invention, the size of the locus L that is selected in each instance of processing 
^'■^ step 204 may be set to a specific size. For example, when the genotypic database 52 

is a SNP database, the size of locus L is set to a predetermined number of 
centiMorgans (cM). Then, in each instance of processing step 204, a different locus 
having the predetermined number of cM is chosen. A centiMorgan is an art 
25 recognized unit of measure that quantifies the spatial relationship between positions 
within a chromosome. More specifically, a centiMorgan is a measure of genetic 
recombination frequency. One cM is equal to a one percent chance that a marker at 
one genetic position will be separated from a marker at another position due to 
crossing over in a single generation. In humans, 1 cM is equivalent, on average, to 1 
30 million base pairs. In some embodiments, the size of the locus L selected in 

processing step 204 is less than 5 cM, 10 cM, 20 cM, 30 cM, 50 cM, 100 cM or a 
value greater than 100 cM. 

It will be appreciated that units other than cM may be used to set the size of 
the locus L selected in each instance of processing step 204. For example, the size of 
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the locus L may be set in units of nucleotides or even kilobases of nucleotides. In one 
embodiment, once the size of the locus has been initially set in a given session, each 
different locus L that is selected in subsequent instances of processing step 204 is 
chosen such that is has the same size as the locus L that was initially selected. 
5 In processing step 206, a genotypic data structure is established for the 

selected locus. In one embodiment, processing step 206 is performed by genotypic 
data structure derivation subroutine 48 (FIG. I). The genotypic data structure is 
typically formed in a method similar to the construction of the phenotypic data 
structure. The values of the phenotypic data structure are typically the differences in 
10 quantitative traits exhibited by several strains of an organism of interest. In contrast, 
the values in the genotypic data structure correspond to counts of the polymorphic 
differences between strains for a given locus L that contains M genetic variations, 
such as SNPs. That is, a given locus L may have several independent genetic 
C3 variations M, and the goal of the genotypic array that corresponds to this locus is to 
U'l 15 quantify the number of these independent genetic variations. To accomplish this, an 
tr individual variation matrix is established for each variation in every position x 
H within locus L. In each such matrix, S"", the i* row and the column are associated 

with the allele value for strain i and the allele value l^'Cj) for strain at locus 
f :^ position X according to the following rule: 

20 

|! S'^O, j) = 1/2 if l^(i) = 0 or m = 0 

= oif r(i) = r(j) 
= iifr(i)7^r(j) 

25 where 0 indicates the allelic value for strain i at locus position x is not known at the 
present time. Therefore, if the alleles for two strains i and j are identical at position x, 
the entry in the individual variation matrix for x would be: 



30 



S''(i,j) = S"(j,i) = 0 

and if the two alleles are different, a "1" is entered. 

In some cases, not all allelic information is known at the present time 
(symbolized by 0). For example, locus position x may contain information on the 
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O 
P 

m 



a 



allele for strain i, but not for strain j. In this situation, the assumption is made that 
strain j has equal probability of containing either allele, and the corresponding entry is 

set equal to one half. 

At this point, in some embodiments of the present invention, each individual 
variation matrix S contains elements that take on one of three values: 0, V2, or 1. It 
will be appreciated that many other types of schemes may be used when allelic 
information is not presently known and use of the value "V2" in such instances merely 
illustrates one example of a scheme that is used in such instances. Similarly, any 
number of weighting schemes can be used rather than a "0" or "1" and all such 
weighting schemes are within the scope of the present invention. 

In one embodiment of the invention, a variation matrix S that ti-acks an 
individual locus position x for five members (Ml through M5) of a species has the 
form: 

Illustrative variation Matrix S 



s 


Ml 


M2 


M3 


M4 


MS 


Ml 


0 


0.5 


0.5 


1 


0 


M2 


0.5 


0 


0.5 


0 


1 


M3 


0.5 


0.5 


0 


1 


1 


M4 


1 


0 


1 


0 


0.5 


M5 


0 


1 


1 


0.5 


0 



In one embodiment of the present invention, in order to assemble the overall 
genotypic matrix for this locus, each individual variation matrix S within the locus L 
selected in processing step 204 is summed. To illustrate this concept, consider the 

20 case in which a locus L was selected in processing step 204 (Fig. 2). In this 

illustrative example, the locus L was selected using a 20 cM window, so the size of 
locus L is 20 cM. Further, there are five locus positions x in locus L. Each locus 
position X is represented by a corresponding variation matrix. In this case, therefore, 
the overall genotypic matrix g(i, j) for this locus is computed by summing the five 

25 variation matrices as follows: 

5 

g(ij) = 2gm(i'j) 



m=5 
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More generally, a given locus L will have M variations, each variation represented by 
a corresponding variation matrix S. Then, the overall genotypic matrix g(i, j) for the 
locus is computed using the expression: 



M 

g(iJ) = ESm(»'j) 

ni=5 



Therefore, an illustrative genotypic matrix G that represents a specific locus in five 
species members (Ml through M5) has the form: 



Illustrative Genotypic Matrix G 



CIS 



G 


Ml 


M2 


M3 


M4 


M5 


Ml 


0 


3.5 


2 


4 


3 


M2 


3.5 


0 


3 


2.5 


1 


M3 


2 


3 


0 


1 


1 


M4 


4 


2.5 


1 


0 


0.5 


M5 


3 


1 


1 


0.5 


0 



10 In viewing the illustrative genotypic matrix G above, it is apparent that there is 

relatively little genotypic variance between members M5 and M4 (0.5) whereas there 
is more variance between Ml and M2 (3.5). 

In one aspect of the invention, each overall genotypic matrix G is assembled 
from individual component variation matrices S within locus L using a weighting 

15 scheme. Generally speaking, a weighting scheme in accordance with the present 

invention first identifies the center of the locus L that was selected in processing step 
204. Variation matrices S that are close to the center of this locus receive full weight 
whereas variation matrices S that are far away from the center of locus L receive only 
partial weight. Thus, the weighting schemes in accordance with the present invention 

20 emphasize or upweight variation matrices S that are near the center of the selected 

locus L and deemphasize or downweight variation matrices that are far away from the 
center of the selected locus L. Weighting schemes in accordance with this aspect of 
the present invention are particularly advantageous when genotypic databases 52 (Fig. 
2) such as SNP databases are used. This is because variation matrices S that are close 

25 to the center of locus L are more reliable than variation matrices S that are far from 
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the center of locus L when such matrices are derived from SNP database data. 
Accordingly, the weighting scheme acts to emphasize more reliable data when the 
data is combined to form genotypic matrix G. 

To illustrate the general principles of the weighting schemes in accordance 
5 with this aspect of the invention, consider the case where a genotypic matrix will be 
generated based upon two variation matrices. Si and S2, that are found within a given 
locus L. 

Si is located 5 cM from the center of locus L and has the values: 



Illustrative variation Matrix Si 



s 


Ml 


M2 


M3 


M4 


M5 


Ml 


0 


0.5 


0.5 


1 


0 


M2 


0.5 


0 


0.5 


0 


1 


M3 


0.5 


0.5 


0 


1 


1 


M4 


1 


0 


1 


0 


0.5 


M5 


0 


1 


1 


0.5 


0 



^ S2 is located 15 cM from the center of locus L and has the values: 

PJ 

\^ Illustrative variation Matrix S2 



s 


Ml 


M2 


M3 


M4 


MS 


Ml 


0 


0.5 


0.5 


1 


0 


M2 


0.5 


0 


0.5 


0 


1 


M3 


0.5 


0.5 


0 


1 


1 


M4 


1 


0 


1 


0 


0.5 


M5 


0 


1 


1 


0.5 


0 



15 

Because S2 is located further away from the center of locus L, one filtering 
scheme in accordance with the present invention applies a weight of 0.5 to each 
element in S2. Therefore, the genotypic matrix G that is derived from the combination 
of all positions x in locus L in this embodiment of the present invention will have the 
20 values: 
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s 


Ml 


M2 


M3 


M4 


MS 


Ml 


0 


0.5 + )^(0.5) 


0.5+%(0.5) 


1+54(1) 


0 


M2 


0.5+ U0.5) 


0 




0 


i+m) 


M3 


0.5+ m.5) 


0.5+>4(0.5) 


0 


i+m) 


1+54(1) 


M4 


1+J^(1) 


0 


1+H(1) 


0 


0.5+ 54(0.5) 


MS 


0 




i+m) 


0.5+>4(0.5) 


0 



It will be appreciated that a broad number of different types of weighting 
schemes may be used to de-emphasize locus positions x that are far away from the 
center of locus L and to emphasize locus positions x that are proximate to the center 
5 of locus L. For example, when genotypic database 52 is a SNP database, the positions 
X in a given locus L can be approximated as a binomial distribution centered on the 
center of locus L. Thus, the distribution of locus positions x about the center of locus 
£3 L may be fitted to a Gaussian probability distribution and each respective locus 
li position X may be weighted by the probability for the respective locus position x that 
S 10 is derived from the Gaussian probability distribution. A Gaussian probability 
p'l distribution weighting scheme is merely provided to demonstrate one form of 
y weighting scheme that is used in some embodiments of the present invention. Many 
t* other forms of weighting schemes based on probability functions are possible. For 
i="^ example Poisson distribution or Lorentzian distribution schemes may be used. See 
P 15 p>>^/^T.atr^n anH Rnhinson. Data reduction and error analysis for the physical sciences, 
McGraw Hill, New York, New York, 1992. 

In some embodiments of the present invention, processing step 206 further 
includes a correlation step in which each gene within the locus L selected in 
processing step 204 is allowed to contribute a maximum of one relative unit to 
20 genotypic matrix G. To illustrate embodiments of the present invention that are in 
accordance with this aspect of the invention, consider the case in which locus L has 
three positions V, d, and V, where D and I? are in gene A and is in gene B. A 
corresponding variation matrix S is computed for each of the three locus positions. 
Then, because each gene is allowed to contribute only one relative unit to the 
25 genotypic matrix G, the variation matrix representing V and the variation matrix 
representing V are given half weight whereas the variation matrix representing d is 
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given full weight when the three variation matrices are summed to yield the 
corresponding genotypic matrix G. 

Embodiments in which each gene within the locus L selected in processing 
step 204 is allowed to contribute a maximum of one relative unit provides an 
5 advantageous filtering effect when correlating phenotypic data to genotypic data in 
subsequent processing steps. Often, in any given genotypic database 52, there are 
some genes that have undergone several mutations and there are some genes that have 
undergone relatively few mutations, if any. After the first few mutations in any given 
gene have arisen, the informative value that subsequent mutations in the gene provide 
10 on localizing phenotypic traits to specific positions in chromosomes diminishes. In 
fact, as the number of mutations in a single gene becomes sufficiently large, the gene 
becomes overrepresented in phenotypic to genotypic correlation computations that are 
performed in the subsequent processing steps illustrated in Fig. 2. To see this, 
consider the case in which a given locus L has two genes A and B and the genotypic 
data for locus L is drawn from a SNP database in which there are ten SNPs for gene A 
and only one for gene B. If genes A and B are not constrained so that they contribute 
one relative unit to the genotypic matrix, gene A would have an order of magnitude 
more influence over gene B in subsequent correlation steps where phenotypic data is 
correlated to genotypic data. This can be seen by an example in which there are two 
strains of mice, Ml and M2, in which genes A and B are represented for Ml and M2 
in a SNP database as follows: 



a 
a 

m 15 

fIJ 

h'^ 20 

o 



Ml: (1, 0); (2, 0); (3, 0); (4, 0); (5, 0); (6, 0); (7, 0); (8, 0); (9, 0); (10, 0); (11, 0) 
M2: (1, 1); (2, 1); (3, 1); (4, 1); (5, 1); (6, 1); (7, 1); (8, 1); (9, 1); (10, 1); (11, 1) 

25 

In the SNP data representation above, each x coordinate represents a position in locus 
L and each y coordinate has a value of "0" when there is a polymorphism present at 
position X and a value of "1" when there is no polymorphism present at position x. In 
this example, positions 1-10 are located in gene A and position 11 is located in gene 
30 B. If genes A and genes B are allowed to contribute unequally to the genotypic 
matrix G, the genotypic matrix will have the values: 
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Genotypic Matrix G 
Unequal contribution from genes A and B 



G 


Ml 


M2 


Ml 


0 


11 


M2 


11 


0 



If genes A and genes B are constrained so that they contribute a maximum amount of 
5 one relative unit to the genotypic matrix, positions 1 through 10 will be weighted by 
0.1 so that they contribute a total of 1. Thus, the genotypic matrix G will have the 

values: 



Genotypic Matrix G 



G 


Ml 


M2 


Ml 


0 


2 


M2 


2 


0 



10 unc 

It 

b 
o 

m 

i^l Imposing the constraint that each gene in a locus L contributes a single relative unit to 

I" the genotypic matrix has the advantage of preventing any given gene or sets of genes 

p j from dominating the correlation coefficient that is computed between phenotypic data 

15 and genotypic data in subsequent processing steps. There are several different ways 
to constrain the relative contribution of each gene within the locus L selected in 
processing step 204 so that any given gene does not overly dominate the 
corresponding genotypic matrix G. For example, genes could be constrained based on 
their length, where longer genes are allowed to contribute more than shorter genes. In 
20 another example, genes could be constrained based on percent A+T nucleotide 

content. In other schemes, genes that include more locus positions x in locus L are 
allowed to contribute more to the genotypic matrix than genes that include less locus 
positions x. However, the amount that such genes are allowed to contribute is not 
linearly proportional to the number of locus positions x within the gene. Rather, for 
25 example, the amount that a particular gene is allowed to contribute to the genotypic 
matrix is logarithmically proportional to the number of locus positions x in the gene. 

In some embodiments of the present invention, two locus positions V and Z/ in 
locus L are considered to be in the same gene if both positions map to a region of 
DNA that has been assigned the same accession number in a genetic database. 
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Genetic databases include databases such as the Human Genome Database (GDB), 
Saccharomyces Genome Database (SGD), Mouse Genome Database (MGD), 
Drosophila Genetic Database (FLYBASE IMGT/LIGM) 

http://www.ebi.ac.uk/embl/DocumentationAJser_manual/dr_line.html) or Genbank 
5 (http://www.ncbi.nlm.nih.gov/Genbank/). Many other genetic databases are known 
and are within the scope of the present invention. 

Now that various embodiments used to construct genotypic matrices have 
been described, attention turns to how these matrices are used. One embodiment of 
genotypic data structure derivation subroutine 48 converts the genotypic matrix to a 
10 genotypic array by taking the non-redundant, non-diagonal elements of the matrix and 
arranging them into the vector G: 

G^giiax gKi,3), giiM. gf<2,3), gK2,4), guM. ...gK/^-i, ^/) 



^5 The vector G obtained for the illustrative genotypic matrix set forth above is G 

g1 = (3.5, 2, 4, 3, 3, 2.5, 1,1,1, 0.5). Once a genotypic matrix such as G has been 

established in processing step 206, a correlation value is formed between the 
phenotypic array and the genotypic array (processing step 208). This correlation 

w 

value is typically computed by phenotype / genotype comparison subroutine 50 (FIG. 
H 20 1). In one embodiment, this correlation is determined by linear regression correlation 
|=a. in which the correlation coefficient is calculated as: 



r(p(i)-<P>) {gii)-<Gf->) 

c(F, G^) = Eqn. 1 

25 {[ 2' - < P >f] [t (00) -<Cf->) 

where, 

c(P, G^) is the correlation value between the phenotypic array and the 
genotypic array that corresponds to locus L; 
30 p(i) is a value of the i'*' element of the phenotypic array; 

g{i) is a value of the i^ element of the genotypic array; 
<P> is a mean value of all elements in the phenotypic array; 
<G^> is a mean value of all elements in the genotypic array; and 
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2' =Z where, 

N is equal to a number of elements in the genotypic array. 
It will be appreciated that the phenotypic and genotypic arrays can be compared in 
processing step 208 using any number of algorithms other than linear regression. For 
5 example, alternative methods for forming a correlation value in processing step 208 
include, but are not limited to, regression analysis, regression analysis with data 
transformations, Pearson correlations, Spearman rank correlation, a regression tree 
and concomitant data reduction, partial least squares, and canonical analysis. {See 
e.g. Lui, "Statistical Genomics," CRC Press LLC, New York, 1998; Stuart & Ord, 
10 "Kendall's Advanced Theory of Statistics," Arnold, London, England, 1994). 

In some embodiments of the present invention, the correlation coefficient is 
weighted by the number of locus positions x in locus L. Such weighting is based on 
p the observation that correlations c(P, G^) computed using a locus L that has a 
i';' relatively large number of locus positions x receive a correlation coefficient that is 
t 15 artificially low relative to those correlations c(P, G^) that are computed using a locus 
5 1 L that has a relatively few number of locus positions x. To illustrate, consider a first 
^ correlation coefficient having the value of 0.5 that was computed using a locus L that 

includes 100 single nucleotide polymorphisms (SNPs) versus a second correlation 
^f '^ coefficient having that value of 0.6 that was computed using a locus L that includes 
P 20 only 10 SNPs. The first correlation coefficient may have more significance because it 
was computed across a much larger number of SNPs. 

It will be appreciated that weighting correlation coefficients c(P, G^) based on 
the number of locus positions x over which they are computed may be performed by 
any number of techniques and all such techniques are within the scope of the present 
25 invention. 

One method of weighting involves computing a correlation coefficient for 
each locus L selected in processing step 204 using the expression: 



30 



[S'(p(i)-<P>) (gKi)-<G^>)][nr 

c(P, G^) = Eqn. 2 

{[ r (p{\) - < P >f] ir (gKi) - < >) 



where. 
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30 



c(P, G^) is the correlation value between the phenotypic array and the 
genotypic array that corresponds to locus L; 

p(i) is a value of the i* element of the phenotypic array; 
g{i) is a value of the i"*" element of the genotypic array; 
<P> is a mean value of all elements in the phenotypic array; 
<gH> is a mean value of all elements in the genotypic array; 
n is the number of locus positions x in locus L; and 

Y =i where. 



i=l 



N is equal to a number of elements in the genotypic array. 
10 It will be appreciated that Eqn. 2 may be derived from Eqn. 1 by multiplying the 

numerator of Eqn. 1 by the square root of n, where n is defined as the number of locus 
positions X in locus L for which a correlation c(P, G^) is being computed. It has been 
determined that, for some data sets, weighting c(F, G^) by the square root of n 
m provides improved c(P, G"-) values. While not intending to be limited to any particular 

0? 15 theory, it is believed that Eqn. 2 corrects for an inherent bias against correlation 
coefficients computed for loci L, using Eqn. 1 , that have a large number of locus 
positions x. Other forms of weighting based on number of locus positions x in locus 
L are possible. For example, rather than multiplying the numerator of Eqn. 1 by the 
square root of n (Eqn. 2), the numerator of Eqn. 1 could be multiplied by n, n^, n 
y, 20 raised to any power. log(n), ln(n), or e°. One of skill in the art will recognize that 

other forms of weighting using n, the number of locus positions x in the locus L, are 
possible and all such weighting schemes are within the scope of the present invention. 
In some embodiments of the present invention, the genotypic database 52 used is a 
SNP database and the number of positions x in the locus L are the number of SNPs in 
25 the SNP database within the given locus L. 

In another embodiment of the present invention, linear regression or weighted 
linear regression is not used to determine a correlation coefficient. Instead, a 
correlative measure cm/ is computed. A correlative measure om/ in accordance 
with this embodiment of the present invention is: 



[r(p(i)-<P>) (g(i)-<G'>)] 
omy(P,d-)= ^^"-3 



{[X'(p(i) -<P>f]}"^ 
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where, 

om^P, G^) is the correlation value between the phenotypic array 
and the genotypic array that corresponds to locus L; 
5 p(i) is a value of the i* element of the phenotypic array; 

g(\) is a value of the i* element of the genotypic array; 
<P> is a mean value of all elements in the phenotypic array; 
and 

is a mean value of all elements in the genotypic array. 
10 While processing steps 202 through 206 have been described with reference to 

linear phenotypic and genotypic arrays, it will be appreciated that the methods of the 
present invention are not limited to the comparison of such arrays. Indeed, any form 
Q of data structure having elements that preserve the information in the above described 



C3 



matrices and arrays may be used. For example, rather than using the genotypic array 



ii1 15 described above, the individual variation matrices can be used. Further, rather than 

^1 using the phenotypic array, a phenotypic distance matrix can be used. 

"'''^ Once a correlation value between the phenotypic data structure and a 

1^ = genotypic data structure that corresponds to a particular locus L has been formed, the 

Til 

1^?^ correlation value is stored in processing step 210 so that it can be subsequently ranked 

20 with the correlation value of each of the other loci that are analyzed. 

P 

Processing step 212 is provided so that the procedure can be repeated in an 
iterative fashion for all suitable loci 54 in genotypic database 52 (FIG 1). Thus, in 
processing step 212, a decision is made whether to test an additional locus by asking 
whether all of the loci present in genotypic database 52 (FIG. 1) have been tested. In 
25 one embodiment, when additional loci 54 are present in genotypic database 52, 
processing step 212 returns a "y^s" and the process continues by looping back to 
processing step 204 where an additional, untested locus is selected from genotypic 
database 52. 

In typical embodiments of the present invention, step 212 acts as a sliding 
30 scale. In such embodiments, an initial instance of processing step 204 picks a locus at 
a starting point on a particular chromosome in the organism of interest. The locus is 
considered a window. This window typically has a length that is measured in 
centiMorgans. Steps 204 through 210 are then performed for the window selected in 



26 



288881.6 



processing step 204. This results in a correlation value for the window. Then, 
process control returns to step 204 where the window is incrementally advanced to a 
position along the chromosome that is contiguous with or even overlaps with a locus 
that was selected in a prior instance of processing step 204. This incremental advance 
may, for example, be a specified number of nucleotides or centiMorgans. When the 
specified number of nucleotides or centiMorgans is less than the window length, it 
follows that successive windows selected in each instance of processing step 204 will 
overlap with each other. The iterative process of selecting a window in processing 
step 204, computing the corresponding correlation value, and advancing the window 
continues until the end of the chromosome is reached. In organisms that have 
multiple chromosomes such as mice, this process continues for each chromosome 
until a window has been advanced over each chromosome in the organism. In one 
embodiment of the present invention, the window is advanced by 10 cM in each 
successive instance of processing step 204. However, this increment is readily 
adjustable. 

In another aspect of the present invention, the window is advanced in each 
successive instance of processing step 204 by a step that approaches an infinitesimally 
small quantity. It has been found that such embodiments provide smoother output. 
Thus, in embodiments where the window is advanced by a very small incremental 
amount, the window is advanced by 2 cM, 1 cM, 0.1 cM, 0.01 cM or less. 

In some embodiments of the present invention, processing step 214 does not 
compute a correlation value using linear regression. Rather, correlative measures 
using equations such as Eqn. 2 or Eqn. 3 are used. The use of a correlative measure 
rather than a correlation coefficient determined by linear regression does not affect 
other aspects of the present invention. 

When there are no additional loci to test (212-No), the correlation value for 
each of the comparisons of genotypic data structures to the phenotypic data structure 
are ranked with respect to each other in processing step 214. In one embodiment, 
processing step 214 comprises the arrangement of the tested loci in a vector K 
according to their correlation scores: 

A:=(LU",r,...) 
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where c(P, G^) > c(P, G^) > c(P, G^O > . . . . 

In another embodiment of the present invention, processing step 214 includes 
the computation of (i) a mean correlation value that represents a mean of each 
correlation value formed during instances of processing step 208; and (ii) a standard 
deviation of the mean correlation value based on each of the correlation values 
formed during instances of processing step 208, 

In processing step 216, the genotypic data structures that achieve the highest 
correlation values are selected. Since each genotypic data structure corresponds to a 
particular locus in the genome, the selection process in processing step 216 results in 
the association of the phenotype with particular loci in the organism of interest. In 
one embodiment, the selection process in processing step 216 is performed by 
selecting genotypic data structures that form a correlation value that is a 
predetermined number of standard deviations above the mean correlation value. 
Typically, the predetermined number is chosen so that a small percentage of the 
genome of the organism, such as five percent, will be selected during processing step 
216. 

In some embodiments of the present invention, phenotype / genotype 
processing module 44 (Fig. 2) includes a user interface. An exemplary user interface 
is illustrated in Figs. 7-10. In some embodiments, the user interface allows the user to 
quickly toggle between a mode in which genotypic matrices are computed in an 
unweighted fashion, where each SNP is given equal weight, and a weighted fashion, 
where each accession number is given equal weight. One of skill in the art will 
appreciate that genotypic data is often characterized by accession numbers, where 
each accession number corresponds to a different gene in an organism of interest. 
Furthermore, in any given genotypic database, there will be several SNPs within any 
given gene. Therefore, each gene, or accession number, will include many SNPs. In 
fact, larger genes will have more SNPs. Thus, weighting by accession number (by 
gene) will produce a very different result then the case where each SNP is given equal 
weight. 

Fig. 7 illustrates a user interface 700 in which a toggle 702 allows the user to 
compute genotypic matrices by accession number. That is, each accession number in 
the locus L selected in processing step 204 (Fig. 2) is given a single "vote" in 
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computing the corresponding genotypic matrix. In Fig. 7 the name of a plurality of 
different mouse strains is listed in panel 704. Further, for each of the mouse strains, 
values for a particular phenotype that correspond to the respective mouse strain is 
shown in panel 706. A panel of check boxes 708 is further provided in user interface 
5 700. The check boxes allow the user to determine which strains will be used in the 
computations of the present invention. Accordingly, when a strain is not selected 
using the check box that corresponds to a given strain, the phenotypic data of that 
strain is not used to compute the phenotypic data structure constructed in processing 
step 202 (Fig. 2). After computations in accordance with Fig. 2 are run, the 
10 correlation coefficient or correlative measure between genotypic data and phenotypic 
data is plotted in panel 710. In panel 710, the x-axis is chromosome location in the 
organism of interest. The y-axis is the number of standard deviations that a particular 
£1 correlation coefficient or correlative measure is above the median correlation 

coefficient or correlative measure from the set of correlation coefficients or 
Ul 15 correlative measures computed using the processing steps disclosed in Fig. 2. For 
instance, peak 712 represents a particular 20 cM window in the genome of a mouse 
- that has a correlation coefficient that is 3.92 standard deviations above the median 

correlation coefficient. Panel 7 10 may be considered a correlation map of the genome 
of the organism under study. 
N° 20 Fig 8. illustrates the same user interface 700 illustrated in Fig. 7. However, in 

Fig. 8, toggle 702 is set so that genotypic matrices are computed by individual SNPs. 
Thus, in the setting shown in Fig. 8, genotypic matrices are computed in an 
unweighted fashion, where each SlSfP gets one "vote" in computing the genotypic 
matrix. 

25 Some embodiments of the present invention provide a user toggle 902 (Fig. 9) 

that allows the user to switch between unweighted and weighted modes. When in 
weighted mode, correlative measures are computed in instances of processing step 
208 (Fig. 2). Each correlative measure is weighted by the number of locus positions x 
within locus L that is represented by the correlative measure. When in unweighted 
30 mode, a correlation coefficient is computed in processing step 208 using an algorithm 
such as linear regression. When in unweighted mode, correlation coefficients 
computed in instances of processing step 208 (Fig. 2) are not weighted by the number 
of locus positions x within the locus L that are represented by the correlation 
coefficient. 
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Some embodiments of the present invention provide a user toggle 1002 (Fig. 
10) that allows the user to set a window size. This window size is used to determine 
the size of the locus L that is selected in successive instances of processing step 204 
(Fig. 2). In one embodiment, window size is measured in centiMorgans. However, it 
5 will be appreciated that other units of measure, such as number of nucleotide bases, 
kilobases, or megabases, are possible. 

Examples 

10 Building a murine SNP database. The methods of the present invention are 

particularly useful in embodiments that make use of genetic information from inbred 
strains of an organism of interest. Thus, a genotypic database 52 was developed that 
fct contains allele information across 15 inbred strains. At Roche Bioscience, 293 SNPs 
O at defined locations were identified in the mouse genome. The SNPs were identified 
Ul 15 by direct sequencing of PCR amplification products from defined chromosomal 
locations. This database also incorporates published allele information for 2848 
'"'^ SNPs, 45% of which are characterized in a subset of M. Musculus strains, and 55% 

of the SNPs are polymorphic between M. castaneus and one or more M. musculus 
f:f subspecies (Lindblad-Toh, et al.. Nature Genetics Apr;24, 381-386, 2000). User 

1==^ 20 queries regarding SNPs found within a specified chromosomal region or between 
selected inbred strains are executed in real time and provided via user interface 24 
(Fig. 1). 

Example 1: Hypothetical example of the method for prediction ofQTL 
regions. To aid in the understanding of the methods of the present invention, FIG. 3 

25 is provided. FIG. 3 shows hypothetical comparisons, in accordance with the methods 
of the present invention, between three mouse strains (A, B, C) using SNP 
information found in the murine SNP database. Each of the two chromosomes sets 
for a given mouse strain is represented by a horizontal box along the horizontal axis 
of FIG. 3. Each chromosome set is characterized by the hatching type (horizontal, 

30 diagonal, and vertical). Chromosomes with the same hatching style in each of the 
mouse strains are identical. Cross hatched or diagonally hatched ovals respectively 
represent alleles at specific chromosomal positions. A dashed horizontal line is used 
to differentiate each of the mouse strains and the accompanying chart at the bottom of 
nG.3. 
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In the hypothetical example provided in FIG. 3, two of the three strains, (A) 
and (B), exhibit a similar phenotype. That is, strains A and B exhibit a similar 
phenotype (full size tail), while strain C has a different phenotype (short tail). SNP 
alleles at particular chromosomal regions are represented as cross hatched or 
5 diagonally hatched ovals. A series of pairwise comparisons, in accordance with the 
algorithm illustrated in FIG. 2, are made to establish the correlation value between the 
phenotype and genotype for each locus. In each of these series of pairwise 
comparisons, allelic differences in a respective segment of the chromosome of each of 
the mouse strains is correlated with the phenotypic difference between each mouse 
10 strain. Graphic analysis of the correlation data between the respective strains is 

shown at the bottom of FIG. 3. The analysis indicates that while most sites exhibit a 
negative correlation with respect to murine tail length, two chromosomal regions 
(302) and (304) have a strong positive correlation. In fact, 302 and 304 are the 
C3 chromosomal regions predicted to have genes regulating tail length. 
U'l 15 The following four examples, (Examples 2 through 5) are made with reference 

^f^ to Fig. 4, Fig. 4 illustrates the correlation between the genotype and phenotype 

^^4 distributions for all 19 mouse autosomal chromosomes for a given trait. Loci are 

arranged proximal to distal for each chromosome. Each bar represents a 30 cM 
^'^ interval of the respective chromosome and neighboring bars are offset by 10 cM. 
20 Dotted line 402 represents a useful cutoff for analyzing the data, with the highest 

fi 

i^l correlated ten percent of the genome being above this line. 

Example 2: Predicting the chromosomal location of the MHC complex. The 
methods of the present invention were used to predict the chromosomal location of 
the MHC complex, which has been mapped to murine chromosome 17, using the H2 

25 haplotypes for the MHC K locus for 10 inbred strains (Anonymous, JAX Notes 475, 
1998). Phenotypic distances for strains that shared a haplotype were set to zero, and a 
distance of one was used for strains of different haplotypes. The SNPs within and 
near the MHC region had a genotypic distribution which was highly correlated with 
the phenotypic distances; the correlation value for interval 440 (FIG. 4A) was 5.35 

30 standard deviations above the average for all loci analyzed. There were no other 

peaks throughout the mouse genome that exhibited a comparable correlation with the 
phenotype. The computational analysis, executed in accordance with the methods of 
the present invention, excluded 96% of the mouse genome from consideration without 
missing the genomic region known to contain the MHC. 
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Example 3: Identification of the QTLs that correspond to allergic asthma. 
The chromosomal positions that regulate susceptibility to experimental allergic 
asthma have been investigated using prior art techniques. For example, pubhshed 
analyses of intercross progeny between susceptible (A/J) and resistant (C3HmeJ) 
5 mouse strains identified QTL intervals on chromosomes 2 and 7 (Ewart, et al., Am J 
Respir Cell Mol Biol 25. 537-545. 2000; Kaip. et al.. Nature Immunology 1, 221-226, 
2000). The ability of the methods of the present invention to identify these 
chromosomal regions was investigated. 

The phenotypic distance used to populate the phenotypic matrix was the 
10 absolute difference between the measured airway response after allergen-challenge 
for each strain pair. The experimentally identified QTL intervals on chromosomes 2 
and 7 were among the strongest peaks identified by the methods of the present 
g invention (FIG. 4B). The computational method excluded 80% of the mouse genome 

from consideration without missing the experimentally mapped QTL regions using 
m 15 airway responsiveness data from only 5 inbred mouse strains. 

Example 4: Lifespan data. Lifespan data for five mouse strains, which 
^■i reflected susceptibility to T cell lymphoma, has been published (Chrisp et al., 
U Veterinary Pathology 33, 735-743. 1996). Using conventional techniques, three 

susceptibility regions were experimentally identified by analysis of intercross progeny 
(Wielowieyski et al.. Mammalian Genome 10, 623-627. 1999; Gilbert, et al., J.Virol. 
67 2083-2090. 1993; Mucenski et al.. Molecular & Cellular Biology 6, 4236-4243. 
1986; Mucenski et al. Molecular & Cellular Biology 8, 301-308, 1988); and all three 
regions were predicted by the computational genome scan (FIG. 4C). In this example, 
over ninety percent of the genome could be excluded from consideration by the 
computational method without overlooking any experimentally verified QTL interval. 

Example 5: Retinal ganglion cells. In another example, the measured density 
of retinal ganglion cells was used as a phenotype. Using conventional techniques, the 
QTLs associated with this phenotype have been localized to chromosome 11 in the 
mouse genome (Williams et al.. Journal of Neuroscience 18, 138-146, 1998). The 
experimentally verified QTL interval on chromosome 1 1 was contained in the 
chromosomal regions predicted by the methods of the present invention, while 96% of 
the mouse genome was excluded (FIG. 4D). 

Example 6: Additional phenotypic traits. The ability of the computational 
method of the present invention to identify candidate chromosomal regions that are 
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associated with six additional quantitative traits was performed. The chromosomal 
positions for these six additional quantitative traits are derived from published studies 
that provided mapped locus positions (quantitative trait loci; QTLs) as well as 
phenotypic data across multiple inbred strains for each trait (Table 2). As shown in 
Table 2, a total of 10 QTLs from 6 published phenotypic studies are identified from 
the literature. Each QTL resides on a different chromosome. Centimorgan positions 
were interpreted from published marker locations on physical maps. 

Table 2: Published chromosomal positions of QTLs that have been associated with 



Phenotype 


Chromosome (cM) 


Notes 


AHR 


2 (23.5), 7(1) 


Allergen induced airway response 
(APTI) 


Eye weight 


5 (0-10) 


Mouse eye weight (grams), day 75 


Retinal anglion 


11 (57.5) 


Retinal ganglion cell # 


Lymphoma 


1 (62-73), 6(30), 16(50) 


Tumor incidence, lifespan 


MHC 


17 (10) 


H2 K serotyping 


PKC 


11 (66), 3(16.4, 45) 


PKC-a protein amount, activity 



The ability of the methods of the present invention to correctly predict 
chromosomal regions containing experimentally verified QTL intervals associated 
with the six phenotypic traits is presented in Table 3. 

Table 3: Summary of predictions made in accordance with 
the methods of the present invention 



Phenotype 


Experimentally 
Verified 


Methods of the Present Invention 


Correct 


Predicted 


Threshold (%) 


AHR 


2 


2 


8 


19 


Eye weight 


1 


1 


6 


17 


Ganglion 


1 


1 


2 


4 


Lymphoma 


3 


3 


4 


8 


MHC 


1 


1 


1 


2 


PKC 


2 


2 


6 


2, 11 


Totals 


10 


10 


27 
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As shown in Table 3, the methods of the present invention identified all ten 
experimentally characterized QTL intervals. In addition, seventeen other 
chromosomal regions were predicted by this computational method. Whether these 
predicted regions affect phenotypic traits has not yet been experimentally verified. 
5 The threshold required for correct identification of a QTL varied from two percent to 
nineteen percent of the complete mouse genome. 

The percentage of correct predictions as a function of the percentage of the 
mouse genome contained within the predicted chromosomal regions was examined. 
If predicted regions contained eighteen percent of the mouse genome (by selecting 
10 eighteen percent of the peaks with the highest correlation), all ten experimentally 

verified QTL intervals were correctly identified (HG. 5). As the threshold was raised, 
limiting the number of predicted candidate chromosomal regions, the methods of the 
present invention missed some experimentally verified QTL intervals for these traits. 
When only three (or nine) percent of the genome was above the threshold, the method 
Si 15 identified four (or seven) of the ten verified QTL intervals for these traits (FIG. 5). 
p When a genome-wide threshold of ten percent was used, the genomic region 

to search for candidate genes was computationally reduced by an order of magnitude. 
Since the average size of a predicted genomic region was 38 cM, the 1500 cM mouse 
genome could be subdivided into approximately forty regions. The computational 
method was used for seven different phenotypes, so approximately 280 genomic 
intervals (38-cM in size) were examined. This method correctly identified seven of 
ten experimentally validated QTL intervals, while missing three, at the ten percent 
genome-wide threshold. The algorithm further predicted 23 genomic intervals were 
involved in a phenotypic trait where no QTL had been experimentally characterized. 
25 Finally, the computational method and experimental analysis agreed on 240 loci that 
were not QTL intervals for the phenotypes examined. This data can be assembled 
into a 2 X 2 matrix to assess the ability of the computational method to predict QTL 
intervals. A Fisher Exact test yields a highly significant P value (7.0 x 10'^) for the 
computationally predicted intervals. 
30 In summary, the methods of the present invention were able to identify ten 

QTLs for seven phenotypic traits that had been previously identified by prior art 
techniques. Each of the experimentally verified QTL intervals was identified by the 
methods of the present invention. The genotypic array used to identify these 
chromosomal regions was derived from a murine SNP genotypic database. In each 
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case, the conventionally identified QTL interval exhibited a computational SNP 
distribution that was highly correlated with the tested phenotype. The correlation was 
well above the mean value for the entire genome, and nine of ten were greater than a 
full standard deviation above the mean. 

5 Example 7: Use of alternative genotypic databases 52. Although the 

examples provided herein utilize a genotypic database of 15 inbred mouse strains, 
other types of genotypic databases may be used. For example, suitable genotypic 
databases include various databases that have various types of gene expression data 
from platform types such as spotted microarray (microarray), high-density 

10 oligonucleotide array (HDA), hybridization filter (filter) and serial analysis of gene 

expression (SAGE) data. 

As a proof of concept, 315 microsatellite polymorphisms were downloaded 

from the Center for Inherited Disease Research URL 
Q http://www.cidr.jhmi.edu/download/CIDR_mouse.xls 

m 15 Genotypic database 52 was populated in manner analogous to the case when SNP data 
^ was used to populate database 52: if the polymorphisms matched between two mouse 

'^J strains, a "0" was entered, if they differed, a "1" was entered. In this way. the number 

U of differences between mouse strains was counted for a given locus. The remainder 

f;l of the analysis was performed in accordance with the methods of the present 

^ 20 invention. For this trial, the MHC locus was identified on chromosome 17. Although 
the QTL for the MHC region was not as clearly distinguished when using 
microsatellite information as it was for SNP data, it should be noted that the 
microsatellite data used for the trial was sparser than the information currently 
available in the mouse SNP database. 
25 Example 8: Comparison of the performance of a genotypic database 52 

populated with SNP data to a genotypic database 52 populated with microsatellite 
data. The genotypic database 52 populated with microsatellite data as described in 
Example 7 was compared to the previously described genotypic database 52 that 
contains allele information across 15 inbred strains for 287 SNPs at defined locations 
30 in the mouse genome. In this case, the phenotype is the formation of retinal ganglion 
cells in infant mice. The experimentally verified QTL that correlates with this 
phenotype is on chromosome 1 1. As illustrated in Fig. 6, the genotypic database 52 
populated with the microsatellite information more strongly identifies the correct QTL 
peak than the genotypic database 52 populated with SNP data (4.2 standard deviations 
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SI 



with microsatellites versus 2.3 standard deviations with SNPs). Furthermore, the 
results using the microsatellite data are less noisy than the results using the SNP data. 
See, for example, the reduced positive peak on chromosome 9 using the microsatellite 

data (602 versus 604). 
5 Example 9. Use of Perturbations. The present invention may be used to 

correlate phenotypes of a plurality of strains of a biological sample with specific 
positions in the genome of the biological sample before and after the biological 
sample has been exposed to a perturbation. In this approach, two sets of experiments 
are performed. In the first set, the methods of the present invention are used to 
10 correlate genotypes to phenotypes before the plurality of strains of the biological 
sample are exposed to a perturbation. In the second set of experiments, the plurality 
of strains of the biological samples are each exposed to a perturbation and the 
p methods of the present invention are used to correlate genotypes to phenotypes. 

Then, the correlations computed in the first set of experiments are compared to the 
U'l 15 correlations computed in the second set of experiments. By comparing differences or 
similarities between these two sets of correlations, it is possible to identify regions of 
^■^ the genome of the biological sample that are highly responsive to the perturbation. In 

U one embodiment of the present invention, the biological sample is a mouse or rat. 

One embodiment of the present invention provides a method of determining a 
1"'^ 20 portion of a genome of an organism that is responsive to a perturbation. In the 
f i method, a first phenotypic data structure that represents a difference in a first 

phenotype between different strains of said organism is produced. The genome of the 
organism includes a plurality of loci. The first phenotype is measured for each of the 
different strains of the organism when each of these different strains is in a first state. 
25 Next, a genotypic data structure is established. The genotypic data structure 

corresponds to a locus selected from the plurality of loci. Further, the genotypic data 
structure represents a variation of at least one component of the locus between 
different strains of the organism. The first phenotypic data structure is compared to 
the genotypic data structure to form a correlation value. These establishing and 
30 comparing steps are repeated for each locus in the plurality of loci, thereby identifying 
a first set of genotypic data structures that form a high correlation value relative to all 
other genotypic data structures that are compared to the first phenotypic data structure 
during the comparing step. 
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The method proceeds with the computation of the second phenotypic data 
structure that represents a difference in a second phenotype between different strains 
of the organism. The second phenotype is measured for each of the different strains 
of the organism when each of the different strain is in a second state. This second 
state is produced by exposing each strain of the organism to a perturbation. 

Next, the second phenotypic data structure is correlated to the genotypic data 
structure to form a correlation value. The computing and correlating steps are 
repeated for each locus in the plurality of loci, thereby identifying a second set of 
genotypic data structures that form a high correlation value relative to all other 
genotypic data structures that are compared to the second phenotypic data structure 
during the correlating step. Finally, a dissimilarity in the first set of genotypic data 
structures and the second set of genotypic structures is resolved, thereby determining 
15 the portion of the genome of the organism that is responsive to the perturbation. 
£^ The phenotypes selected for study in the two sets of experiments may be any 

m 15 type of phenotype that is reliably measured. Thus, a phenotype may be, for example, 
life-span of the biological sample, the basal serum level of an antibody in the blood of 
y the biological sample, the serum level of an antibody in the blood of the biological 
U sample after exposure of the biological sample to a perturbation, the response of a 

biological sample in one of the various pain models described in Example 10 after the 
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20 biological sample has been exposed to a pain relieving drug, etc. Many other 



25 



30 



phenotypes are possible and all such phenotypes are within the scope of the present 
invention. 

The term "perturbation" within the context of this example is broad. A 
perturbation can be the exposure of a biological sample to a chemical compound such 
as a pharmacological or carcinogenic agent, the addition of an exogenous gene into 
the genome, or the removal of an exogenous gene. Thus, for example, the antibody 
serum level in mice representing a plurality of difference mice species can be 
measured before and after exposing each strain of mice to an antigen. Then, the 
genotypic differences in the plurality of different mouse strains is correlated with 
observed phenotypes before and after exposure of the mice to a perturbation. By 
comparing the peaks found in the correlation map of the mice before and after 
exposure to the perturbation, it is possible to localize regions of the mouse genome 
that are most affected by the perturbation. 
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The panel of check boxes 708, provided in user interface 700 (Fig. 7), is 
particularly useful in cases where perturbations are used. For any given perturbation, 
there will typically be a strain having a phenotype that is more responsive to the 
perturbation than all other strains studied. To determine how the highly responsive 
strain affects the correlation map plotted in panel 7 10 of Fig. 7, one simply deselects 
the unresponsive species and reruns the calculations. 

Once the regions of the genome that are highly responsive to the perturbation 
have been identified, gene chip expression libraries that include the identified portion 
of the genome may be examined. Of particular interest is the identification of 
differential expression of genes in (i) a gene chip library made from a strain of the 
biological sample before insult with a perturbation and (ii) a gene chip library made 
from the strain of the biological sample after insult with a perturbation. As is well 
known in the art, the gene chip library may be a collection of mRNA expression 
levels or some other metric, such as protein expression levels of individual genes 
within the organism. Comparison of the differential expression level of genes in the 
U two gene chip libraries leads to the identification of individual genes that exhibit a 
5 high degree of differential expression before and after exposure of the biological 
\ sample to a perturbation. Correlation of the positions of these individual genes with 

K the regions of the genome identified using the correlation metrics disclosed above 

ll 20 provides a method of identifying specific genes that are highly responsive to a 
perturbation. 

Exemplary gene chip expression libraries have been used in studies such as 
those disclosed in Karp et al. "Identification of complement factor 5 as a 
susceptibility locus for experimental allergic asthma," Nature Immunology 1(3), 221- 
226 (2000) and Rozzo et al. "Evidence for an Interferon-inducible Gene, Ifi202, in the 
Susceptibility of Systemic Lupus," Immunity 15, 435-443 (2001). Furthermore, 
methods for making several different types of gene chip libraries are provided by 
vendors such as Hyseq (Sunnyvale California) and Affymax (Palo Alto, California). 

Example 10. The following protocols illustrate some of the many ways in 
which phenotypic data can be derived for biological samples of interest, in order to 
practice the methods of the present invention. 
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1. IN VIVO ACTIVITY IN RATS^ 

The following protocols are generally as described in Faden, 1989, Brain 
Research 486:228-235 and Mcintosh et al., 1989, Neuroscience 28(l):233-244. 

1.1 Animals used. Male Sprague-Dawley rats (375-425g) are obtained from 

5 Harlan (Frederick, MD) and housed for at least 1 week prior to any procedures. The 
animals are maintained at a constant temperature (22 ± 2°C) and a 12 hr light/dark 
cycle, with lights on at 6 am and all neurological scoring performed during the light 
cycle. Food and water are available ad libitum. 

1.2 Fluid-Percussion Induced Traumatic Brain Injury (TBI). Rats are 

10 anesthetized with sodium pentobarbital (70 mg/kg i.p.), intubated, and implanted with 
femoral venous and arterial catheters. Brain temperature are assessed indirectly 
through a thermister in the temporalis muscle. Body temperature is maintained 
1^ through a feedback-controlled heating blanket. Blood pressure is continuously 

P monitored, and arterial blood gases analyzed periodically. After the animal is placed 

m 15 in a stereotaxic frame, the scalp and temporal muscle are reflected, and a small 
^ craniotomy (5 mm) located midway between the lambda and bregma sutures over the 

0 

\l left parietal cortex allows insertion of a Leur-Loc that is cemented in place. The 

fluid-percussion head injury device, manufactured by the Medical College of 
I Virginia, consists of a plexiglass cylindrical reservoir filled with isotonic saline; one 

1,1 20 end includes a transducer that is mounted and connected to a 5 mm tube that attaches 
fj through a male Leur-Loc fitting to the female Leur-Loc cemented at the time of 

surgery. A pendulum strikes a piston at the opposite end of the device, producing a 
pressure pulse of approximately 22 msec duration, leading to deformation of 
underlying brain. The degree of injury is related to the pressure pulse, expressed in 
25 atmospheres (atm): 2.6 atm in our laboratory produces a moderate injury with regard 
to neurological and histological deficit. Sham (control) animals undergo anesthesia 
and surgery without fluid percussion brain injury. 

L3 Neurological Scoring. Standardized motor scoring is performed at 1, 7 
and 14 days after TBI, by individuals unaware of treatment. Motor function is 
30 evaluated utilizing three separate tests, each of which is scored via an ordinal scale 
ranging from 0=severely impaired to 5=normal function. Tests include ability to 
maintain position on an inclined plane in the vertical and two horizontal positions for 
5 sec; forelimb flexion (suspension by the tail) and forced lateral pulsion. Each of 
seven individual scores (vertical angle, right and left horizontal angle, right and left 
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forelimb flexion, right and left lateral pulsion) are added to yield a composite 
neurological score ranging from 0 to 35. This scoring method shows high interrater 
reliability and is very sensitive to pharmacological manipulations (see, Faden et al., 
1989, Science 244:798-800). 

5 1.4 Automatic and Analeptic Assessment. Additional groups of uninjured rats 

are tested for autonomic and analeptic responses immediately prior to and up to 60 
minutes following drug administration. For the analeptic study, rats are first 
anesthetized with 40 mg/kg i.p. sodium pentabarbitone and placed onto an unheated 
pad on the laboratory bechtop at room temperature (22 ± 2°C). A thermister probe is 

10 placed in the rectum to measure core body temperature. After a ten minute period, 
rats are administered vehicle or drug as described below via the tail vein. Time to 
recovery of the righting reflex was subsequently determined while temperature is 
recorded at five minute intervals for all animals. 

To assess autonomic responses to perturbations, such as pain relieving drugs, a 

15 separate group of rats are anesthesized with 4% isoflurane (1.5 L/min). Catheters are 
then placed into the right artoid artery and right jugular vein and exteriorized at the 
back of the neck. Rats are separated one per cage and allowed to recover from 
anesthesia. The exteriorized catheters are suspended above the rat to prevent biting. 
Mean arteriolar blood pressure (MAP) is continuously recorded via a transducer 

20 connected directly to the arterial catheter for the duration of the study. At 1 h 

following catheter placement, each rat is administered vehicle or drug via the catheter 
in the jugular vein as described below. 

1.5 Administration of Compounds. Rats are injected via the femoral vein 
catheter with a single bolus dose (1 mg/kg) with various compounds of interest. The 

25 investigator is blinded to drug treatment both at the time of surgery and for 

neurological scoring. For autonomic and analytic studies, rats are given either normal 
saline or a compound under study at the times indicated above. 

1.6 Data Analysis, Continuous variables compared across groups are 
examined using an analysis of variance (ANOVA) followed by Bonferroni correction 

30 (rightin reflex). Continuous variables subjected to repeated measurements over a 
period of time (cardiovascular and core temperature measurements) are analyzed 
using a repeated measurements ANOVA followed by Tukey's pairwise comparison at 
each time point. Ordinal measurements (composite neurological scores) are evaluated 
using the non-parametric Kruskal-Wallis ANOVA with individual, non-parametric 
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Mann-Whitney U-tests. Survival differences are compared using the Chi-Square test. 
A p value <0.05 is considered statistically significant. 

2. IN VIVO STUDIES IN MICE 

5 2.1 Animals. Male C57B1/6 mice (20-25g) are obtained from Taconic Farms 

(Germantown, NY) and housed in an area directly adjoining surgical and behavioral 
rooms for at least 1 week prior to any procedures. All mice are maintained at a 
constant temperature (22±2°C) and a 12 hr light/dark cycle, with lights on at 6 am and 
all behavioral testing performed during the light cycle. Food and water are available 

10 ad libitum. 

2.2 Controlled Cortical Impact Device. The injury device consists of a 
microprocessor-controlled pneumatic impactor with a 3.5 nmi diameter tip. The 
impactor is vertically mounted on a mill table (Sheriine, USA) which allows for 
precise adjustment in the vertical plane above the mouse head, which itself is secured 

15 to a stereotaxic apparatus (David Kopf Instruments, CA) attached to the instrument. 
The core rod of a linear voltage differential transducer (LVDT, Serotec, USA) is 
attached to the lower end of the impactor to allow measurement of velocities between 
3.0 and 9.0 m/s. Velocity of the impactor is controlled by fine tuning both positive 
and negative (back) air pressures. An oscilloscope (Tektronix, USA) records the 

20 time/displacement curve produced by the downward force on the LVDT, allowing 
precise measurement of the impactor velocity. 

2.3 Surgery, Surgical anesthesia is induced and maintained with 4% and 2% 
isoflurane respectively, using a flow rate of 1.0 - 1.5 1 oxygen per minute. Depth of 
anesthesia is assessed by monitoring respiration rate and palpebral and pedal- 

25 withdrawal reflexes. The animal is then placed onto a heated pad and core body 

temperature is monitored and maintained at 38 +/- 0.2 "^C. The head is mounted in a 
stereotaxic frame and the surgical site clipped and prepared with a series of three 
Nolvasan scrubs followed by sterile saline rinses. A 10 mm mid-line incision is made 
over the skull, the skin and fascia reflected, and a 4 mm craniotomy made on the 

30 central aspect of the left parietal bone with a tissue punch (Roboz, USA). Great care 
is taken with the removal of the parietal bone to avoid injury to the underlying dura 
mater which is continuously bathed in sterile normal saline warmed to 37.5 °C. The 
impounder tip of the pneumatic injury device is cleaned with a pad, soaked in 
absolute alcohol, positioned to the surface of the exposed dura and automatically 
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withdrawn the 44 mm stroke distance. Following injury at a moderate (6.0 m/s 
velocity, 1 mm tissue deformation depth) level, the incision is closed with interrupted 
6-0 silk sutures, anaesthesia is discontinued and the mouse was placed into a heated 
cage to maintain normothermia for 45 minutes post-injury. All animals are monitored 
5 carefully for at least 4 hours post-surgery and then daily. To minimize variation 
between animals due to anaesthesia during acute neurological testing, 20 minutes is 
allowed for surgery and five minutes for suturing for each animal. 

2.4 Administration of Compounds. Conscious mice are placed in a mouse 
restrainer and injected via the lateral tail vein with either normal saline or a compound 

10 of interest at 30 minutes following controlled cortical impact injury (CCI). The 
investigator is blinded to drug treatment both at the time of surgery and for 
neurological and behavioral scoring. 

2.5 Acute and Chronic Neurological Evaluation. Chronic neurological 
1^1 recovery is evaluated for all animals using a beam walking task, a method that is 

15 particularly good at discriminating fine motor coordination differences between 

injured and sham-operated animals. The device consists of a narrow wooden beam 6 
? j mm wide and 120 mm in length that is suspended 300 mm above a 60 mm-thick foam 

^ rubber pad. The mouse is placed on one end of the beam and the number of footfaults 

11 j for the right hindlimb recorded over 50 steps counted in either direction on the beam. 

20 A basal level of competence at this task was established before surgery with an 
O acceptance level of <10 faults per 50 steps. 

2.6 Spatial Learning Evaluation. The Morris watermaze (Morris, 1984, J. 
Neurosci. Meth. 22:47-60) is employed to assess spatial learning by training mice to 
locate a hidden, submerged platform using extramaze visual information. The 

25 apparatus consists of a large, white circular pool (900 mm diameter, 500 mm high, 

water temperature 24 ± 1°C) with a plexiglass platform 76 mm diameter painted white 
and submerged 15 mm below the surface of water (225 mm high) that is rendered 
opaque with the addition of dilute, white, non-toxic paint. During training, the 
platform is hidden in one quadrant 14 cm from the side wall. The mouse is gently 

30 placed into the water facing the wall at one of four randomly-chosen locations 

separated by 90 degrees. The latency to find the hidden platform within a 90 second 
criterion time is recorded by a blinded observer. On the first trial, mice failing to find 
the platform within 90 seconds are assisted to the platform. Animals are allowed to 
remain on the platform for 15 seconds on the first trial and 10 seconds on all 
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subsequent trials. There is an inter-trial interval of 30 minutes, during which time the 
mice are towel-dried and placed under a heat lamp. A series of 16 training trials 
administered in blocks of 4 are typically conducted on days 7, 8, 9, and 10 post- 
surgery. 

5 Example 11. Effect of constraining each gene to a single vote. The advantage 

of constraining each gene in a locus L to a single vote when constructing a genotypic 
matrix from component variation matrices S will now be disclosed. Fig. 1 1 shows a 
correlation map where each variation in each locus selected in successive instances of 
processing step 204 is allowed to contribute to the corresponding genotypic matrices 
10 irrespective of whether multiple variations exists in a single gene. Thus, in the 

computation of the correlation map 1 102 of Fig. 11, multiple SNPs in the same gene 
contribute to the corresponding genotypic matrix if they fall within the locus selected 



O in processing step 204 (Fig. 2). The data in panel 1 102 is a plot of the correlation 
ft coefficient computed between respective genotypic and phenotypic arrays across the 
^Is entire mouse genome. The correlation map shows a peak 1 104 that is 2.8 standard 
h deviations above the mean correlation score computed for the entire map. The gene 
J that is known to affect the trait under study in Fig. 1 1 is actually in chromosome 17 at 

15 cM. Thus, the peak in Fig 1 1 is in the wrong region of the mouse genome. 
|.a In Fig. 12, each gene in the locus L selected in processing step 204 is 

f 320 constrained to a single vote in a parliamentary style. Thus, if there are multiple 
variations in the particular gene, each variation is scaled so that the sum of the 
variations equals a single vote. When this form of constraint is imposed, the 
correlation map across the entire mouse genome reveals a peak 1202 that is centered 
on a gene that is known to influence the trait under study. Furthermore, the peak is 
25 now 4.05 standard deviations above the mean score. 



Discussion 

Computational analysis of genotypic databases 54 using phenotypic data from 
sources such as inbred parental strains and the methods of the present invention 
30 rapidly identifies candidate QTL intervals. This can eliminate many months to years 
of laboratory work required for generation, characterization and genotyping of 
intercross progeny. In effect, the methods of the present invention reduce the time 
required for QTL interval identification from many months to milliseconds. 
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There are several factors contributing to the successful QTL predictions by 
computational scanning of the murine SNP genotypic database using the methods of 
the present invention. The use of inbred mouse strains limits variability due to 
environment, and timed experimental intervention and sampling limits error in 
5 phenotypic assessment. The inbred strains are homozygous at all loci, which 

eliminates confounding effects due to heterozygosity found in human populations. 
However, there is no absolute requirement that inbred strains be used to populate 
genotypic database 52. 

The methods of the present invention will greatly accelerate analysis of 
10 complex traits and mammalian disease biology. Recently, there has been increased 
emphasis on using chemical mutagenesis in the mouse as a method for studying 
complex biology. This has occurred as a result of the difficulties noted by 
investigators searching for complex trait loci using standard methods for QTL 
Ff analysis. For a review, see Nadeau and Frankel, Nature Genetics Aug;25, 381-384, 
W is 2000. However, analysis of genetic variation among existing inbred mouse strains 
can be markedly accelerated by application of the methods of the present invention. 
""'^ Of course, understanding the genetic basis of complex disease requires additional 

steps beyond computational prediction of genomic intervals. Specific gene candidates 
must be identified and evaluated before the underlying mutations can be identified 
^'^20 and effective treatment strategies can be designed, tested in animal models, and 

o 

developed for use with humans. 

Alternative Embodiments 

The foregoing descriptions of specific embodiments of the present invention 
25 are presented for purposes of illustration and description. They are not intended to be 
exhaustive or to limit the invention to the precise forms disclosed, obviously many 
modifications and variations are possible in view of the above teachings. For 
example, the techniques of the invention may be applied using pooled or clustered 
genetic variation information as a source for the genotypic data structure or genetic 
30 variation information from individual samples. Similarly, the phenotypic information 
provided from sources, such as phenotypic data file 60, may be in the form of pooled 
or clustered phenotypic data or phenotypic data from individual organisms. 
Furthermore, genotypic database 52 may represent inbred strains of the organism of 
interest or randomized strains of the organism of interest that have not been inbred. 
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4 • 

Because of the overwhelming homology between murine and human genomes, the 
examples provided herein clearly demonstrate that the methods of the present 
invention provide an invaluable tool for correlating human phenotypic traits with 
specific loci in the human genome. 

While the examples provided herein describe the comparison of a plurality of 
genotypic data structures to a phenotypic data structure, one of skill in the art will 
appreciate that many other types of comparisons may be practiced in accordance with 
the present invention. For instance, consider the genotypic to phenotypic data 
structure comparison as a two-dimensional comparison. Higher dimensional 
comparisons than the two-dimensional comparison are possible. For instance, one 
embodiment of the present invention provides for a three dimensional comparison of 
the class: "genotypic data structure" versus "phenotypic data structure one" versus 
"phenotypic data structure two." Another example of a type of comparison within the 
scope of the present invention includes a comparison of "SNP genotypic data" to 
"disease phenotypic data" to "microarray data." 

References cited and conclusion 

All references cited herein are incorporated herein by reference in their 
entirety and for all purposes to the same extent as if each individual publication or 
patent or patent application was specifically and individually indicated to be 
incorporated by reference in its entirety for all purposes. 

The embodiments were chosen and described in order to best explain the 
principles of the invention and its practical application, to thereby enable others 
skilled in the art to best utilize the invention and various embodiments with various 
modifications as are suited to the particular use contemplated. It is intended that the 
scope of the invention be defined by the following claims and their equivalents. 



45 



288881.6 



